Cyclotron resonance lineshape in a Wigner crystal 
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The cyclotron resonance absorption spectrum in a Wigner 
crystal is calculated. Effects of spin-splitting are modelled 
by substitutional disorder, and calculated in the coherent po- 
tential approximation. Due to the increasing strength of the 
dipole-dipole interaction, the results show a crossover from a 
double-peak spectrum at small filling factors to a single-peak 
spectrum at filling factors > 1/6. Radiation damping and 
magnetophonon scattering can also influence the cyclotron 
resonance. The results are in very good agreement with ex- 
periments. 

PACS numbers: 78.20.Ls, 78.66.-w, 73.20.Dx 



Two-dimensional (2D) electron systems in strong mag- 
netic fields have been studied intensively for the last 15 
years. Most of the work has been concentrated on trans- 
port properties. Important information can also be ex- 
tracted from measurements of high-frequency properties 
such as photoluminescence and the cyclotron resonance. 
A large number of experimental and -theoretical studies 
have dealt with these phenomenal!! □ The interactions 
with the low-energy degrees of freedom influence the de- 
tailed structure of the high-energy resonance. 

Experiments on 2D electron systems at small Landau 
level (LL) .-filling factors v by Besson et alB and Sum- 
mers et ala revealed intriguing changes in the cyclotron 
resonance absorption spectra as the filling factor v, and 
temperature were varied. For very small v the spectra 
show two peaks. With increasing v, the absorption peak 
at the lower frequency increases in relative strength, and 
the peaks are shifted. Eventually, for v > 1/6 the two 
peaks merge. The spectra also change in qualitatively 
different ways with increasing temperature for v smaller 
and larger than m 1/10, respectively. It was speculated 
that this signaled a phase transition, for example the for- 
mation of a Wigner crystal. n 

In a recent paper, Cooper and Chalkertl introduced a 
model that explains the experimental results without in- 
voking a phase transition. The idea behind the model 
is that, due to band-structure effects, spin-up and spin- 
down electrons have slightly different cyclotron resonance 
frequencies.B At a finite temperature there are electrons 
with both spins, and the 2D system is a disordered mix- 
ture of the two spin species. In this model, the elec- 
trons (i.e. their guiding centers) occupy the lattice sites 
of a Wigner crystal (WC). The optical properties of this 
system near the cyclotron resonance frequency can be 
described in terms of excitations, "excitons", from the 
lowest LL to the next LL. The excitons propagate on a 



triangular lattice with substitutional disorder. The prop- 
agation is caused by the electromagnetic dipole-dipole 
interaction, and its strength relative to the splitting be- 
tween the spin-up and spin-down cyclotron frequencies 
determines the shape of the absorption spectrum. In Rcf. 
pL the averaging over different disorder configurations was 
done numerically. 

This paper presents a study of the effects of spin- 
splitting disorder treated within the coherent potential 
approximation (CPA),El which makes it possible to ob- 
tain results that agree with experiment, with only a small 
calculational effort. We also present a thorough discus- 
sion of the physics behind the results. Moreover, we 
deal with two other processes that broaden the cyclotron 
resonance: radiation damping and thermal motion of 
the electron guidin^centers treated in terms of magne- 
tophonons (MP's)BEl Radiation damping gives a contri- 
bution to the resonance width comparable to experimen- 
tal values once the electron density is large enough. In the 
case of a double-peak spectrum, the MP's broaden the 
high-frequency peak more than the low-frequency peak. 

The model that we will study describes excitons prop- 
agating in a WC with substitutional disorder. The WC is 
a convenient model system, but the results of the calcu- 
lations should be similar also in an electron liquid. The 
Hamiltonian can be written 
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where Cj annihilates an exciton at site i. The site cy- 
clotron frequency u>i equals the bare cyclotron frequency 
lu c = eB/m* (B is the magnetic field and m* the electron 
effective mass) fax-Xhe majority spins, and u> c + 5lo for the 
minority spins .qEj The exciton- hopping matrix elements 
due to the dipole-dipole interaction are given by 
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As will be seen, this fixes the 



cyclotron frequency at eB/m* in a WC without disor- 
der, in accordance with Kohn's theorem.til In Eq. (0), 
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are lattice vectors of the WC, eo is the dielectric 



constant of vacuum, e r is the relative dielectric constant 
of the material (GaAs), and l c — (h/eB) 1 / 2 is the mag- 
netic length. The Fourier transform t q of the hopping 
matrix elements yields the exciton dispersion in a system 
without disorder. It is proportional to the trace of the 
dynamical matrix-of a Wigner crystal in the absence of 
a magnetic fieldj3'EI and can be related to the transverse 
and longitudinal phonon frequencies UTq and wt q by 



&(c4_ + wi_)/(2w ) 



Lq) 



(3) 



For large magnetic fields and small filling factors Eq. (g) 
is, for all practical purposes, equivalent to the magneto- 
plasmon dispersion relation resulting from a lattice dy- 
namics treatment. a 

The absorption probability of a circularly polarized 
photon with frequency u> near the cyclotron frequency, 
impinging on the electron system at right angle, is 



P{w) = 



e n e uj c 
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This result treats the A ra( j ■ j interaction between the ra- 
diation, and the electrons to lowest order (Fermi golden 
rule) £3 In Eq. (H), n e is the electron density, cq is the 
speed of light in vacuum, and G(q, u>) is the Fourier trans- 
form of the retarded exciton Green's function 



Gy(t) = -W(t)<[c i (t),ct(0)]). 



(5) 



We use the CPA to calculate G(q,w). The CPA is 
known-,to give correct results in a-^aumbcr of important 
limits|3 and Persson and Rybergc3 used it successfully 
to solve the formally equivalent problem of optical ab- 
sorption in an isotopically mixed adsorbate layer. The 
disordered system, where the exciton frequency at a site 
iSji^c with probability (1 — c) and uj c + 8oj with probability 
c,Ej is replaced by a translationally invariant system (ef- 
fective medium) where the exciton frequency is uj c + (j(uj) 
everywhere. The exciton self energy cr(u)), which must be 
determined self-consistently, is frequency-dependent and, 
in general, complex. The Green's function can be written 



G(q, uS) = [u + «7 — u> c — cr(u)) — tq/h] 



(6) 



where 7 represents exciton damping due to other mech- 
anisms than disorder (see below). In the CPA, one real 
site (with frequency lo c or lo c + Slu) is placed in the ef- 
fective medium of io c + o~(uo) sites. Requiring that the 
scattering off the real site, treated to all orders, should 
vanish on the average yields an equation for a(ui) J3 



(j{u) 



cSlo 



1 - [§u> - er(w)]GooM 
The diagonal Green's function in real space 
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where Abz is the area of the first Brillouin zone (BZ). In 
the calculations, we have evaluated Eq. (@) accurately, 
but the analytic approximationMI 
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where a = lu + ij — lo c — cr{ui) and W is the exciton 
bandwidth, gives almost the same results. 



There are a number of processes that can cause ex- 
citon damping described by the constant 7 in Eq. (jg). 
Radiative decay is one of the most, important, and it can 
be calculated in a standard way.Ej Excitons with a wave 
vector larger than ^Je^oj c /cQ do not radiate. For smaller 
wave vectors the radiation damping is 
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(neglecting exciton dispersion). This result resembles 
those for the decay of free excitons in quantum wells.E^I 
Numerically 7i7o=0.013 meV when n e = 10 11 cm~ 2 . 

It should be emphasized that Eq. (|J) yields P < 1 with 
radiation damping included, whereas neglecting it may 
give P > 1. On resonance u> = w c , and without disorder 
crjui) = 0, setting 7 = 7 q in Eq. (g), and then using Eq. 
(0) yields P = 1. The present calculations suppress the 
q dependence of the damping. The constant 7 = 7 q =o + 
7 P h, using hj p h = 0.005 meV for the phenomenological 
damping constant. 

Figure [j] shows the main results of this calculation. In 
Fig. [y(a) the magnetic field and temperature are kept 
constant, while the electron density varies. The quali- 
tative behavior of the spectra as a function of n e is in 
close agreement with experiment (cf. Fig. 2 in Ref. [|). 
The bare cyclotron resonance lies at hcu c — 27.64 meV 
and the minority spin resonance is at 27.49 meV. If n e 
is chosen small enough, the absorption spectrum has two 
isolated resonances at these frequencies. However, all the 
spectra displayed in Fig. |l](a) are modified by the dipole- 
dipole interaction, since it is comparable in strength to 

the difference in resonance frequency between the two 

3/2 
spin species. The exciton bandwidth W ~ n e increases 

from 0.10 meV at v = 1/14.2 to 0.67 meV for v = 1/4.1. 

The physics behind the results of Fig. 0(a) can be un- 
derstood in the following way (cf. Refs. [RJ and |5|): The 
transition between the two lowest LL's at a site can be 
considered as a local oscillator. The high-frequency, ma- 
jority oscillators respond in phase with the incident exter- 
nal radiation field when it has a frequency near the lower 
resonance. Then, the total field acting on the minority 
oscillators is enhanced due to anti-screening, and the low- 
frequency peak grows in strength with increasing n e . At 
the same time, it is shifted towards higher frequencies be- 
cause the majority-oscillator response is larger there. At 
the high-frequency peak, on the other hand, the major- 
ity oscillators absorb energy from the external field and 
thus have a phase lag of w tt/2. The minority oscillators 
now respond out of phase with the total field and there- 
fore lead the external field. Consequently, the minority 
oscillators effectively emit radiation, and this reduces the 
strength of the second peak in the absorption spectrum. 

Through all the changes of the absorption spectrum 
its first moment stays constant at uj m = lo c + c5lo. In 
the end, at high densities, only one peak centered at ui m 



remains. This behavior is characteristic of motional nar- 
rowing. Once the hopping matrix elements become large 
enough, exciton hopping becomes so frequent that the 
excitons are not much affected by the disorder. In this 
regime, the disorder contribution to the linewidth be- 
haves as ~ Suj 2 /W, and decreases with increasing den- 
sity. In an impurity-scattering language, the factor 1/W 
comes from the density of final states. The other contri- 
butions to the linewidth are not influenced by this nar- 
rowing, and in Fig. [y(a) the linewidth starts to grow at 
the highest densities due to radiation damping. This re- 
sult is also in agreement with experiment .Q 

Figure pl(b) shows absorption spectra at different tem- 
peratures and minority spin concentrations. At low 
enough temperature, there are only majority spins, and 
consequently only one peak. When the temperature is 
increased, this peak splits in qualitatively different ways 
depending on the density. 

Consider now the effects of electron thermal and zero- 
point motion. The calculation of the exciton-hopping 
matrix elements in Eq. (0) assumes that the electrons are 
fixed at the lattice sites. In reality they perform vibra- 
tions so that the matrix elements tij are time-dependent. 
To quantify this time-dependence, return to Eq. (0) and 
replace the lattice site coordinates Rj(j) by Rj(j) +u i ( J ), 
where Ui(j) is the electron ^displacement that can be de- 
scribed in terms of MP'sB Then expand the resulting 
expression to linear order in Uitj) recovering the static 
exciton-hopping matrix element plus a new term, describ- 
ing exciton-MP interaction due to anharmonicity. This 
exciton-MP interaction Hamiltonian can be written 
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where a q annihilates an MP of wave vector q. 

The expression determining the matrix element M qq ' 
is lengthy however, M qq / = whenever q or q' or their 
difference vanishes. Thus, in a system without disorder, 
as a direct consequence of Kohn's theorem, a q = ex- 
citon is not scattered by the MP's. But the exciton-MP 
interaction still influences the cyclotron resonance line- 
shape in the presence of disorder. The disorder scatters 
a long-wavelength exciton into an intermediate state with 
a finite q, from where the exciton can decay due to the 
interaction with the MP's. 

We include the effects of these processes on the exci- 
ton spectrum through a two-step calculation. First, the 
decay rate of an exciton with wave vector q, caused by 
H cx _ mp , is calculated in a system without disorder. The 
Fermi golden rule gives 

r,= yE l M <N'! 2 ft 1 + ^q-q')*(*q - *q' - &"q-q') 
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(ph , and the same calculations as before yield the exciton 
self energy and Green's function. 

Figure g shows the results of these calculations. The 
most interesting behavior occurs for an intermediate 
value of the electron density (n e — 3.2 x 10 10 cm" 2 ). The 
low-frequency peak is slightly broadened when MP scat- 
tering is included, but this is a minute change. The high- 
frequency peak, on the other hand, is substantially broad- 
ened, and only a weak shoulder remains. This seems to 
be in qualitative agreement with the experimental results 
of Ref. 0. There (see Fig. 1), for a similar filling factor, 
the high-frequency peak is more temperature-sensitive 
than the one at lower frequency. Similar results were also 
found in Ref. |[ where these thermal effects were mod- 
elled by frozen structural disorder. At lower densities, 
when the spectrum approaches that of two isolated reso- 
nances, both peaks are somewhat broadened by the MP 
scattering. No results for the motional narrowing regime 
are shown, since the MP scattering has little influence on 
the spectrum there. 

The above results can be explained by very much the 
same reasoning as was used earlier. In the intermedi- 
ate density regime and at frequencies near the first peak, 
the oscillators respond more or less in phase with each 
other. The small changes in the dipole-dipole interaction 
caused by the MP's do not then perturb the spectrum ap- 
preciably. For higher frequencies, near the second peak, 
nearest-neighbor oscillators may be out of phase with 
each other. Changes in the coupling due to the MP's 
are much more important in this case. 

In conclusion, this paper has presented a model calcu- 
lation of the cyclotron resonance lineshape in a Wigner 
crystal. The effects of spin-splitting of the cyclotron res- 
onance has been treated within the CPA. The result- 
ing absorption spectrum has two relatively independent 
peaks, one for each spin species, at low electron densi- 
ties. When the density increases, more and more strength 
is transferred to the low-frequency peak due to anti- 
screening effects until, for filling factors v > 1/6, only 
one peak remains. The radiation damping contribution 
to the linewidth was calculated and found comparable to 
experimental linewidths for v > 1/4. Finally, a study of 
the effects of magnetophonon scattering showed that it 
has small effects on the spectrum, except under special 
circumstances. All of the results are in good agreement 
with experiments. 

I have benefited from interesting discussions with Bo 
Persson on several occasions. I thank Alan Luther for 
useful comments on the manuscript. 



where w q is an MP frequency and JV q is the thermal MP 
occupation. Then ij is replaced by 1(7 + r q /2) in Eq. 
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FIG. 1. Calculated absorption spectra, (a) The electron 
density is varied as indicated next to the curves, while the 
magnetic field and temperature are kept constant. The minor- 
ity spin concentration is 10 %. At the lowest density, the two 
spin species give one resonance peak each, but anti-screening 
effects transfer spectral strength from the high-frequency, to 
the low-frequency peak. Eventually, with increasing density, 
one motionally narrowed peak emerges. It starts to broaden 
at the highest densities because of radiation damping, (b) 
Spectra for two different electron densities n e (3.8xl0 10 cm -2 
and 2.7xl0 10 cm -2 ), with varying temperature (and minor- 
ity spin concentration) as indicated next to the curves. The 
spectra split in qualitatively different ways depending on n e . 



FIG. 2. Absorption spectra calculated with and without 
exciton-magnetophonon (MP) scattering. The MP's are not 
effective in scattering long-wavelength excitons. Thus, for 
the higher electron density, only the high-frequency peak is 
considerably affected by the MP's. At the lower density, both 
peaks are slightly broadened. 
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